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The low-frequency shot noise properties of two electrostatically coupled semiconductor quantum 
dot states which are connected to emitter/collector contacts are studied. A master equation ap- 
proach is used to analyze the bias voltage dependence of the Fano factor as a measure of temporal 
correlations in tunneling current caused by Pauli's exclusion principle and the Coulomb interaction. 
In particular, the influence of the Coulomb interaction on the shot noise behavior is discussed in 
detail and predictions for future experiments will be given. Furthermore, we propose a mechanism 
for negative differential conductance and investigate the related super-Poissonian shot noise. 
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I. INTRODUCTION 

Shot noise investigations in mesoscopic systems can re- 
veal information of transport properties which are not ac- 
cessible by conductance measurements alone [jj. In par- 
ticular, the dynamic correlations in the tunneling current 
through double-barrier structures caused by Pauli's ex- 
clusion principle can provide information regarding the 
barrier geometry 0. If the charging energy of bound 
states becomes larger than the thermal energy as in the 
case of small quantum dots (QDs), strong Coulomb cor- 
relations occur and have an additional influence on the 
shot noise. For metallic QDs the zero-frequency Fano 
factor which quantifies correlations with respect to the 
uncorrelated Poissonian noise |3( was analyzed by a mas- 
ter equation approach including the Coulomb blockade 
effect in Ref. jij- At the steps of the resulting Coulomb 
staircase the Fano factor shows dips caused by Coulomb 
correlations which is quantitatively confirmed in the ex- 
periment In Ref. |(| a similar theoretical approach 
was applied to determine the finite-frequency shot noise 
of metallic QDs. 

Despite that, for semiconductor QDs a comprehensive 
picture of the bias dependent shot noise behavior and 
a subsequent comparison with experimental data is not 
available at the moment. Some theoretical work has 
been done: the investigation of the Fano factor of an 
ensemble of states with statistically varying positions in 
a barrier by a classical approach pj; analytical bias de- 
pendence of the Fano factor for zero-temperature with 
non-equilibrium Green's functions yah bias dependence of 
spin-dependent coherent tunneling [9j; shot noise in the 
co-tunneling regime |1C|: in the Kondo- regime In 
Ref. [12| a numerical investigation of a spatially extended 
QD by means of a coherent technique was presented. In 
the current plateau regime the authors find a suppressed 
Fano factor which is smaller than one half for symmetric 
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barriers in contradiction to the result of |2(. In contrast, 
an enhanced Fano factor at the current steps was found. 
All of the above references mainly consider the noise due 
to negative correlations (sub-Poissonian noise). In the 
case of negative differential conductance in the current- 
voltage characteristic, e.g. in resonant tunneling diodes, 
positive correlations lead to super-Poissonian noise [lflj . 
Furthermore, in capacitively coupled metallic QDs super- 
Poissonian noise can occur |14| . 

Recently, the measurement of the low-frequency shot 
noise of tunneling through an ensemble of self-organized 
QDs was presented in Ref. [l5| which primarily motivates 
this work. The current- voltage characteristic for low bias 
is dominated by steps which are presumably due to tun- 
neling through few QD ground states which are well sep- 
arated in energy (see also HHHGl). 

The correspond- 
ing Fano factor shows an average noise suppression on 
the current plateau which enables the determination of 
the effective collector barrier thickness, given the known 
thickness of the emitter barrier At the current steps, 
Fano factor peaks appear which we have considered the- 
oretically by a master equation approach 0. It was 
shown that for tunneling through QD states which are 
not subject to Coulomb interaction, Fano factor peaks 
at the bias position of current steps occur, caused by 
Pauli's exclusion principle. Good qualitative agreement 
with experiment was found. 

The goal of this paper is the exemplary demonstration 
of the influence of the Coulomb interaction upon the 
Fano factor in a system of two QD states, its interplay 
with Pauli's exclusion principle, and the consequences for 
future noise experiments where Coulomb interaction is 
present. The outline of the paper is as follows: in Se c. UTI 
a brief description of the master equation formalism [2(j 
and the calculation of spectral power density (where we 
mainly follow the lines of Ref. (3) is given. Sec. IHII con- 
tains the results of the bias dependent Fano factor for 
varying Coulomb interaction energy and in Sec. IIVI the 
super-Poissonian shot noise related to negative differen- 
tial conductance will be discussed. All results will be 
summarized in Sec. [V] 
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ensemble average. A certain state v = (m, . . . ,njv) (oc- 
cupation numbers rii G {0, 1}) at a time i of the con- 
sidered QD system with N single-particle states is de- 
scribed by the respective occupation probability P v (t). 
The single-particle states may correspond to different en- 
ergy levels in the same QD, or in different QDs. The time 
evolution of these occupation probabilities is determined 
by sequential tunneling of an electron into or out of the 
emitter/collector contact with the tunneling rates T E j c 
(i labels the single-particle QD state) and can be written 
as a master equation in the general form 



Figure 1: Scheme of the tunneling structure with rates ^e/c 
(i = 1,2) between two QDs and the emitter/collector, respec- 
tively; single-particle levels Ei/%, and Coulomb interaction 
energy U 



P = MP 



(2) 



The spectral power density of the current fluctuations is 
related to the autocorrelation function of the current by 
the Wiener-Khintchinc theorem: 



S o6 H = 2 / dte™* ((I a (t)I b (0)) ~ (if) U) 



where a,b = E,C (emitter/collector) and (•) denotes the 



Note that this approach holds only for weakly coupled 
QDs, 

^e/c ^ k E T, and cannot account for co-tunneling 
processes, i.e. coherent tunneling of two electrons simul- 
tanously [ToL l2*l| . In this paper we consider two QD 
states Ex 1 2 m two different QDs which are connected 
to the emitter/collector contact and are coupled elec- 
trostatically by Coulomb interaction of strength U (see 
Fig. Then the vector of the occupation probabilities 

is £ = ( P (o,o)> ^(1,0^(0,1) !-P(i,i)) T and the transition 
matrix in J2J) reads [2(| 
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T E and the Fermi functions 



with r := T l c + + 

in the emitter p E = (l + exp((J5; - er]V) / '(ksT)))- 1 and 
ff = (l + exp((Ei + U -eijF)/^)))- 1 . V is the 
bias voltage, r]V is the voltage drop across the emitter 
barrier, and f E U includes the Coulomb interaction en- 
ergy U of the occupied QDs. For eV 3> k E T we neglect 
tunneling from the collector into the QDs, setting the 
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The steady state solution of © can be obtained by 
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so that the stationary mean current reads 
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In order to calculate the current flowing through the sys- 
tem current operators are introduced. The current oper- 
ators at the collector barrier with p c = f^, u = 
the emitter barrier, respectively, are defined by 



fp, U = and at 



In the stationary limit the current at the collector bar- 
rier equals the mean current at the emitter barrier. For 
the calculation of the stationary current in JfJJl the cur- 
rent operators JSJ could also be defined in diagonal form 
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without changing the result (JSJ. But for the determi- 
nation of the current-current correlator (see below) the 
definition J^J becomes crucial as it projects the occupa- 
tion probability to the state after an electron traversed 
the barrier. 

To define the autocorrelation function of the current the 
time propagator T_(t) is introduced as follows: 

T(i) =; exp (M*) with P(i) =T(t)P(0) (7) 

With I0J , (0 , and the current-current correlator in 
is 3 



k B T = AE/23 and 7 =5). The mean current (/) vs. 
bias voltage V is plotted in Fig. [SJi. For U = there are 
two steps due to tunneling through the respective states. 
The width of the current steps is determined by the Fermi 
distribution of the emitter electrons. Note that typical 
energy scales are as follows: the bias voltage V is of the 
order of tens of mV, AE can be of the order of a few 
meV, kpT is of the order of tens of fieV for temperatures 
in the range of a few Kelvins. 

The respective Fano factor a © is shown in Fig. [2b. On 
the first plateau in the current-voltage curve of Fig. [2b 
where only tunneling through one single-particle state 
occurs the Fano factor becomes 



(I a (t)I b (Q)) = *(*)£[] T(t)j P°]„ + 

£ — ' =a — =b 



+eS ab 8(t) £ 



[j P°]„ 

—ajb 



(8) 



(0(t) is the Heaviside function). The first two terms of 
the right-hand side in (JSJ contain the correlation between 
tunneling events at different times. The last term de- 
scribes the self-correlation of a tunneling event at the 
same barrier (for further discussions of |JHJ see 0,H2,|2j|). 
For uj <C (* = 1>^) wn i cn i s the regime where 

experimental data are available at the moment (e.g. 
|15|). the spectral power densities become constant and 
See(0) = S CC (0) = S EC (0) = S CE (0) = S(0) holds. 
As a measure of deviation from the uncorrelated Poisso- 
nian noise the dimensionless Fano factor a is used Q: 



a(0) = 



S(0) 



Sp 

with S P = 2e(I) 



(9) 



(Poissonian noise 



III. COULOMB INTERACTING QUANTUM 
DOTS 

We consider the tunneling through two non-degenerate 
QD ground states E\ and E2 with an energy separation 
AE = E2 — E\ which could correspond to slightly differ- 
ent QD sizes. The coupling to the emitter and collector 
contact will be assumed to be the same for both states: 

r R = - ^ 



and r c 
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In the following we discuss three cases for the Coulomb 
interaction energy U: A. U = 0, B. U < AE, and C. 
U > AE. 



A. Noninteracting states: (7 = 

In Figs. |21 and |3 the results of a calculation for variation 
of U in the range of a few ksT are shown (for fixed 
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with 



where i = 1 denotes the energetically lowest state. For 
f E — 1 eq. I|1U|I is the well-known relation derived by 
L. Y. Chen et al. 2]. It reflects the sensitivity of the 
Fano factor to Pauli's exclusion principle. This is shown 
by the full curve a(j) in Fig. For symmetric tunnel- 
ing barriers a is equal to one half and approaches unity 
for strong asymmetry. For bias voltages below the cur- 
rent onset where f\ w the tunneling current becomes 
uncorrelated so that a.\ = 1. At the second step where 
the second single-particle state is filled the Fano factor 
(Fig. |2b) has a peak which is also an effect of Pauli's 
exclusion principle and we obtain a simple analytical ex- 
pression for an arbitrary number of noninteracting QD 
states (for a derivation for two states see Appendix A): 



(I) 



(11) 



by (fTU|) . where the current through state i is: (/j) = 

ej-^-f E , and the net current is (I) = Y^iiW- Eq- l(TT|) 
was applied to the measured Fano factor modulation 
of tunneling through self-organized QDs 15] in a bias 
regime where only a few QD ground states are active in 
transport. It can qualitatively reproduce the measured 
Fano factor dependence upon the bias voltage 0] . 

B. U < AE 

With increasing Coulomb interaction (7^0 the Fano fac- 
tor peak vanishes (see Fig. [2b) while the current changes 
only slightly. This underlines again the strong sensitivity 
of shot noise to correlations. 

Further increase of U leads to an additional step whose 
bias voltage is proportional to U. The respective occupa- 



tion probabilities Pnj 
are shown in Fig. [zj;: 



0)i °(0,l)j 



and P (1)1) for U = 16k B T 
at the first plateau the electrons 
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Figure 2: a) Normalized mean current (7) vs. bias voltage V . 
b) Fano factor a vs. bias voltage V for different values of the 
Coulomb interaction energy U. c) Occupation probabilities 
^(1,0); ^(0,1)! an d P(i,i) vs - bias voltage V for U/ksT =16. 
Parameters: k B T = AS/23, 7 = V C /T E =5. 



tunnel through the energetically lowest state (1,0); the 
second plateau is generated by tunneling through both 
single-particle states with different probabilities and with 
lower probability through the two-particle state which is 
determined by the coupling to the collector. This cor- 
related state originates from aligning the emitter Fermi 
energy with the energy E\ + AE of the second single- 
particle state which can be filled then. If the system is 
in the state (0, 1), a second electron may enter the i = 1 
level, as U < AE. In contrast, the level i = 2 is not acces- 
sible from the state (1, 0) as long as rjV < E\ + AE + U . 
This explains the asymmetry between the occupation 
probabilities P(i,o) an d -P(o,i)- The height of this sec- 
ond plateau depends on the ratio of the tunneling rates: 
oc (1 + with 7 := T C /T E - 

At the second plateau in the current- voltage characteris- 



tic the Fano factor differs from the case of U = 0, where 
only Pauli's exclusion principle plays a role. The depen- 
dence on the ratio of tunneling rates 7 is shown by the 
dotted curve in Fig. 21 In contrast to the noninteract- 
ing regime the minimum is now at 7 = 2 with a = 5/9. 
These are numerical values, as we did not obtain an an- 
alytical expression for this Coulomb correlated state. 




(1,0) (0,1) 




(eriV-E^/AE 

Figure 3: Same as Fig. |5] but for different values of U. c) 
U/k B T =37. 

In the experiment of Ref . [l5( the question arises whether 
the QD states which are contributing to transport are 
Coulomb interacting. One way of determining this ques- 
tion would be the analysis of the Fano factor dependence 
upon the tunneling rate ratio 7 as shown in Fig.0J How- 
ever, in the experimental setup of Ref. 0| these rates 
are determined by the growth procedure. Therefore, they 
cannot be varied in the same sample. Here, we propose 
how to obtain the information about Coulomb correla- 
tions via the temperature dependence of the Fano factor. 
In Fig.0the Fano factor a vs. bias voltage V for different 
temperatures T in the bias range of the second current 
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Figure 4: Fano factor a vs. ratio of tunneling rates 7 = 
Tc/Te- Full curve: first current plateau in Fig. Dashed 
curve: second current plateau; Dotted curve: third current 
plateau. U = 37fc s T. 



step of Fig. [2] is plotted. For nonintcracting QD states 
(Fig. Oi) the Fano factor peak gets broader and experi- 
ences a slight shift to lower bias voltages for increasing 
temperatures. A qualitatively different picture results for 
interacting QD states in Fig. 03 (U = AE/46): with in- 
creasing temperature the peak increases and also shifts to 
lower voltages. Hence, a unique fingerprint of Coulomb 
interaction of QD states even for very small U shows up 
in the temperature dependence of the Fano factor peaks. 
Experimental investigations of this are in progress |24j . 

C. U > AE 

For U > AE a fourth step arises in the current vs. bias 
voltage characteristic in Fig. Ek). Now, the second cur- 
rent plateau corresponds to a different state as in section 
B. Due to U > AE only the two single-particle states 
can be filled with the same probability (compare Fig- EJ;) . 
The respective Fano factor dependence upon 7 is (as dis- 
cussed in Ref. 0) 

a = 1 - -= — E with f E = 2T E (12) 

and is shown by the dashed curve in Fig. 21 The effect of 
Coulomb correlation upon the Fano factor consists sub- 
stituting by 2T E in eq. H10J1 and leads to a shift of the 
minimum of the full curve in Fig. 01 by 7 = 2. For 7 =1 
the Fano factor is a =5/9. 

IV. SUPER-POISSONIAN NOISE 

In the previous section we used identical rates Tc/e for 
both states. Now, we allow their couplings to the collec- 
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Figure 5: Temperature dependence of the Fano factor vs. bias 
voltage V for 7 =5. a) U = 0. b) U = AE/46 



tor contact to be different: ^ Tq. For 72 < 0.457i 
this leads to negative differential conductance (NDC) in 
the current-voltage characteristic as shown in Fig. [fJJi 
(AE = 0) 25] . If the tunneling rate to the collector 
of state i — 2 is much smaller than the other one the oc- 
cupation probability -P(o,i) becomes close to unity. The 
current carried by this state is proportional to Tq and 
therefore low. Consequently, the occupation probability 
Ph,0) is low except close to the current onset where the 
respective current peak arises and around the bias volt- 
age where the two-particle state (1, 1) becomes occupied, 
(compare Fig. [(J). Let us consider the dependence of the 
Fano factor on the bias voltage in Fig.^Js: in the bias volt- 
age range where the current is suppressed the Fano factor 
is larger than unity (super-Poissonian shot noise caused 
by positive correlations of tunneling events). This is sim- 
ilar to the situation in a resonant tunneling diode [T^ |. 
The Fano factor exhibits two different values separated 
by a step in the middle of the NDC-region (dashed curve 
in Figs.EJ}) and c)) marked by arrows. The Fano factor 
dependence on 72 in the regime where P(i,o) decreases or 
keeps constant is shown in Fig.[7| This is the actual NDC 
regime: for 72 — > a approaches 3. In the bias regime 
where -P(i,o) starts to increase due to thermally activated 
electrons (arrows in Fig. EJ:) the Fano factor is lowered 
due to the effect of Pauli's exclusion principle. Its value 
depends on the coupling ratio 71 . 
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Figure 6: a) Normalized mean current (7) vs. bias voltage V . 
b) Fano factor a vs. bias voltage V. c) Occupation probability 
P(i,q) vs - bias voltage (logarithmic scale) for 72 = (dash- 
dotted curve), 72 = 10 -6 (dashed curve), 72 = 0.001 (full 
curve). Parameters: 71 = 1, k B T = Z7/23, AE =0. Note 
that all three curves coincide in (a). 
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Figure 7: Fano factor a vs. 72 at (erjV — Ei)/U =0.2 in 
Fig. Eb- 



Coulomb correlated in tunneling experiments through 
self-organized QDs. This can be done by investigating 
the temperature dependence of the Fano factor vs. bias 
voltage at steps in the current-voltage characteristic. 
Furthermore, we have examined an NDC mechanism in a 
system where the two states are coupled to the collector 
with different tunneling rates. Then the weakly coupled 
state blocks the other state by the Coulomb interaction. 
For degenerate states a current peak occurs and we have 
analyzed the Fano factor which becomes larger than unity 
(super-Poissonian noise due to positive correlations) in 
the bias range where the current is Coulomb blocked. 
In the limit of vanishing coupling of one state to the 
collector, we find a =3. 
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V. CONCLUSIONS 

We have investigated the low-frequency shot noise behav- 
ior in tunneling through two non-degenerate QD states 
which interact electrostatically. For noninteracting states 
the respective non-equilibrium current- voltage character- 
istic shows steps due to resonant tunneling through the 
single-particle states. In this case we have derived an 
explicit analytical expression for the bias dependence of 
the Fano factor. For interacting states additional current 
steps occur which are associated with Coulomb correlated 
states. The influence of the couplings to the collector 
upon the Fano factor has been clarified. We have de- 
veloped sensitive tools to determine whether states are 



Appendix A: ANALYTICAL EVALUATION OF 
THE FANO FACTOR FOR TWO 
NONINTERACTING STATES 

The stationary probability that level i is occupied is pi = 
^f|a or unoccupied 1 - Pl (T l :=r l E + r l c ). Then the 
stationary occupation probability of the noninteracting 
two-level system given by (@J reads 



(1- Pl )(l-P2) 
-P2) 
P2O- ^-Pl) 
PlP2 



(Al) 
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since in the uncorrelated case the occupation probability 
for each state factorizes into the occupation probabilities 
of the single levels. By inserting this vector into JJjJ one 
immediately sees that terms with PiPj cancel and the 
current is the sum of the currents through each level i: 
(Ii) = eT l c pi. This also holds for an arbitrary number of 
levels: </} = £iW- 

Now, let us consider the time propagator Q: its matrix 
elements T„ M (£) describe the conditional probability to 
have state v at time t under the condition of state \x at 
t = 0. The matrix element T ufl (t) = T^ u {t) can be 
factorized for each level i with the following conditional 
probabilities: 

n, = 0->l: Pl (l-e- rlt ) 

th = l-»0: (1 - Pi )(l - e~ TH ) 

m = 0^0: l-pi(l -e- TH ) 

rii = Pi + e- rH {l- Pi ) (A2) 



Due to the form of the current operator at the collector 
barrier in JSJ the first row and last column of the matrix 
T(t) does not enter in the calculation of the current- 
current correlator JSJ. Carrying out the sum in iJSJ for 
two levels leads to 



(I c (t)I c (0)) = 2eT 1 c {[(h) Pl (l-p 2 ) + (I 2 )p 2 (l -pi)][r (0 ,o)^( 1 ,o)+r (0 ,o)^( 1 ,i)] 

+{h)pi[T(ifi)->(i,o) + T{\, o)— (i,i)] + (-fi)P2p(o,i)->(i,o) +^(o,i)->(i,i)]} 
+ 2eIo -pa) + (^2)^2(1 -Pi)][r (0 , )^(o,i) + 7 1 (o,o)^(i,i)] 

+ (J r 2)Plp 1 (l, 0)^(0,1) + T(l, 0)-»(l,l)] + (^l)P2[T'(0, 1)-»(Q,1) +^(0,1)^(1,1)]} 

+ e(I)S(t) 



Replacing the in (|A3|) by using the rules (|A2|) the 

correlator becomes 

(/o(t)/c(0)> = - 2(/ 1 ) 2 e- rl * - 2(I 2 ) 2 e- r2t + (I) 2 

+ e(I)S(t) (A4) 

which can be generalized for an arbitrary number of levels 
(Ic(t)I c (P))= -2 5>} 2 e - rIt + </) 2 

i 

+ e{I)5(t) (A5) 



(A3) 

I 

The time-independent term in (|A4|) and (|A5|) cancels out 
in the calculation of the spectral power density JJJ and 
we obtain 



S(0)=2e(/)-4^^f (A6) 

i 

Dividing Eq. I|A6(I by 2e(J) and using the Fano factor for 

tunneling through a single level i, i.e. a, — 1 — ^pr- (TDJ, 
the Fano factor for an arbitrary number of noninteracting 
levels eq. (fTT|) is derived. 
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